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Abstract 

We have derived the dipolar relaxation function for a cluster model whose 
volume distribution was obtained from the generalized maximum Tsallis 
nonextensive entropy principle. The power law exponents of the relaxation 
function are simply related to a global fractal parameter a and for large time 
to the entropy nonextensivity parameter q. For intermediate times the relax- 
ation follows a stretched exponential behavior. The asymptotic power law 
behaviors both in the time and the frequency domains coincide with those 
of the Weron generalized dielectric function derived from an extension of the 
Levy central limit theorem. They are in full agreement with the Jonscher 
universality principle. Moreover our model gives a physical interpretation of 
the mathematical parameters of the Weron stochastic theory and opens new 
paths to understand the ubiquity of self-similarity and power laws in the re- 
laxation of large classes of materials in terms of their fractal and nonextensive 
properties. 

1 Introduction 

These last decades, growth, aggregation and fragmentation models have been 
proposed and discussed with the aim of representing the (multi)-fractal ge- 
ometry and the resulting scaling properties of a great variety of complex 
physical systems such as glasses, polymers, colloids, gels, self-similar porous 
and cellular materials P] . These studies cover a large range of physical (elec- 
trical, dielectric, optical, magnetic and mechanical) properties j2]. It has 
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also been observed that universality is the major feature of their scaling 
behavior. 

In the same context, it is more and more accepted that the Boltzmann- 
Gibbs statistics is not adequate to describe the macroscopic thermodynamic 
properties of natural phenomena when the effective microscopic interactions 
and the microscopic memory are long ranged due to complex non-equilibrium 
growth, aggregation or fragmentation processes. For that reason a general- 
ized form of a nonextensive entropy known as the "Tsallis entropy" [S] has 
been used with some success although the fundamental basis of such a formu- 
lation is still the object of many discussions. Recently, Abe and Rajagopal [Hj 
have shown that Gibbs theorem for canonical ensemble theory is not univer- 
sal and that counting rule for statistical description can depend on physical 
systems. This nonuniqueness can lead to various forms of the entropy and to 
appropriate modifications of the maximum entropy principle. An increasing 
number of papers are dealing with this problem [7j jH] . 

Sotolongo et al. [HI starting from Tsallis'expression have used the princi- 
ple of maximum entropy with some natural constraints to derive a fragment 
size distribution function in a regime in which long-range correlations be- 
tween all parts and scaling are present. This very general principle gives rise 
to scaling without additional assumptions and the resulting size distribution 
can be consider as a paradigm of physical systems where multiscale interac- 
tions play the key role in determining the hierarchical form of the fragment or 
cluster size distribution. It is the reason why we have found interesting to use 
this distribution to analyze the relaxation properties of complex materials. 

Relaxation (dielectric, mechanical, magnetic.) in systems such as the 
ones we are considering here is a complex stochastic mechanism which de- 
pends both of the cluster geometric structure and the collective nature of 
the interactions. The whole process of regression towards a steady pseudo- 
equilibrium state when an external field is applied to or release from the 
material is dynamical in nature. In most cluster models, relaxation is viewed 
as a complex hierarchical or stochastic process [TU]|llj . 

We will view here a cluster as a set of units (atoms, molecules, aggregates) 
relaxing collectively due to their interactions and their geometric structure. 
The scale of the relaxation is time dependent. Small clusters relax first as 
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a result of intra-cluster interactions. For larger time due to long-range in- 
teraction relaxing clusters are connected to the already relaxed clusters and 
therefore depend on intra- and inte- cluster interactions. 

We will consider the relaxation process from a macroscopic perspective 
and introduced two global parameters one defining the time and space frac- 
tal nature of the relaxation (a) and a second derived from the maximization 
of the nonextensive entropy (g) and characterizing the hierarchical structure 
of the cluster geometry. In that way we aim at deriving a macroscopic "uni- 
versal" relaxation function valid for many physical systems with different 
relaxation mechanisms at the micro and meso scale. We consider here only 
electrical dipolar relaxation. 

Indeed in that case, an universal pattern independent of materials and 
microscopic model has been noticed for a long time jH],|l] and widely dis- 
cussed. Experimentally one observes both for large t and small t a fractional 
power law dielectric response in the time domain : 

^^'^ = -| = 

[ {upt)-""-^ for ujpt»l 
with < n,m < 1. Quite generally it is observed that 1 — n < m. The fre- 
quency Up is material dependent. For intermediate t, the relaxation function 
(f){t) is usually fitted to a stretched exponential, the so-called Williams- Watts 
form : (f){t) ~ exp{ — {ujpt)^ . The small t (large u) exponent n is related to 
small clusters relaxation and the large t (small frequency) exponent m — 1 
describes the relaxation of large clusters. As it will appear in this work, 
they are not independent. This non-Debye behavior has its counterpart in 
the frequency domain where generalizations of the Cole and Cole dielectric 
functions with exponents linked to n and m are commonly used by experi- 
mentalists. 

The purpose of this letter is to show that our (a, q) model leads naturally 
to this universal behavior and relates simply the exponents n and m to the 
two parameters a and q. Asymptotically and numerically our results coincide 
with the general dipolar relaxation function derived by Weron and collab- 
orators ^T] using generalizations of the Levy central-limit theorem. This 
unexpected result gives a physical meaning to the mathematical parameters 
appearing in the purely stochastic Weron theory. 
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2 Cluster size distribution function 



Following references PIjCII we start from the expression of the generalized 
entropy proposed by Tsallis [S] and formulate the maximization process as 
in 

Sq = k,'-^Jl^^ (2) 

The integral runs over all admissible values of the quantity x and x + dx 
and p{x)dx is the probability of the system being in the state between x and 
X + dx. The quantity g is a real number characterizing the nonadditivity 
of the entropy (nonextensivity) for system having nonexponential statistical 
distribution and fc^ is the Bolzmann constant. When q ^ 1, one recovers 
the BG statistics i.e Sg S = —ks p{x) \np{x)dx. 

If we denote the volume of a cluster by V and some typical volume Vm 
characteristic of the distribution (for example the volume having the linear 
dimension of a disorder or correlation characteristic length), we can define 
a dimensionless volume v = V/Vm which corresponds to the quantity x in 
eq.(2). The first constraint to impose is the normalization condition 

p{v)dv = 1 (3) 





One has then to impose a mass conservation. As the system is finite and 
in order to yield a slow decay in the asymptotic behavior of the cluster size 
distribution, a more general condition, the "q-conservation of the mass" has 
been imposed in the form 

J^jp'^{v)dv ^ ^ 
J^pi{v)dv 

and re- worked on the basis of the "normalized q-expectation values" intro- 
duced in jUj US] . 

Eq.(3) and eq.(4) are the constraints to impose in order to derive the 
g-dependent CSDF (cluster size distribution function) using the method of 
Lagrange multipliers . The extremization calculation leads to the following 



4 



hierarchical CSDF defined for 1 < g < 2: 



p{v) = [1 




2-q J 



(5) 



which has the same formal form as in [01,^2] but with different coefficients 
since here we have used the more correct normalized g-expectation values 
using the so-called "escort" probabilities jlS]. This distribution has been used 
successfully to account for the transition to scaling observed in the behavior 
of fragments in the process of breaking [Hj as well as more subtle effects as the 
dimension cross-over between small and large fragment regions when thick 
clay plates and glass are fractured [T2j. This is the cluster distribution we 
will use to discuss dipolar relaxation in self-similar structured systems. 

3 Cluster relaxation function 

We will assume following a common practice that the observed macroscopic 
non-exponential relaxation behavior can be interpreted in terms of a weighed 
average of a Debye exponential relaxation decay exp[— t/r] with respect to a 
distribution of relaxation time Pt{t)- In this formulation the relaxation func- 
tion that describes the polarization decay with time after a steady polarizing 
electrical field has suddenly been removed, is written as: 



To relate the cluster distribution eq. (5) to the relaxation rate distribution, 
we have to make some assumptions on the variation of the relaxation time 
with the number of relaxing units N within a cluster. In a previous paper on 
relaxation in percolation clusters fHl, following a model of linear polymers, 
we assumed, as a first approximation, that the cluster relaxation time r was 
proportional to the number of atoms (or the volume) in .the percolating 
cluster and obtained in that case the proper long time scaling behavior. 
For more complex systems due to intra- and inter-cluster interactions and 
reasons discussed in the introduction, it is reasonable to assume that the 
relaxing time increases faster than linearly with A^. The relation of r with 
N has been mainly discussed for the stress relaxation in the phenomena of 
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polymerization and sol-gel transition. For example, in a simple Rouse model 
|17j the cluster relaxation time has been identified with the time scale for 
diffusion of the cluster over a distance corresponding to its own size. Since 
in such model = R%/Dn wheieDN ~ and Rn is the radius of 

gyration which scales with the fractal mass of the cluster as Rjy ~ N^^^f , 

we get the relation r = N'^^'^f . Although each relaxation mechanism has its 
own specificity and is material dependent, it is legitimate to assume quite 
generally that the relaxation time of a volume made of N relaxing elements 
scales as: 

r = t;^/" (7) 

with < a < 1. The exponent a represents macroscopically the "fractal" 
geometric and dynamical nature of the relaxation dynamics in other words 
the geometrical and dynamical "exploration" to use an expression coined by 
de Gennes^^ We will therefore assume a relaxation time distribution of the 
form: 

P.M = [1-^^1^ (8) 
z q 

To make the comparison with the Weron JT] stochastic theory easier we will 
write the relaxation function in terms of the distribution of relaxation rate 
(3 = 1/r. We therefore write the relaxation function 0Q,^g(t)after the proper 
change of variable as a Laplace transform relation: 

POO 

4>a,g{t) = / ]9/3(a,g)exp[-t/3]d/? (9) 
Jo 

with 

pp{a,q) = - (10) 

A q 

This normalized function belongs to the basin of attraction of the Levy dis- 
tribution with a Levy tail exponent a. Weron et aZ. have argued that 
the P in eq.(lO) are not individual rates but effective relaxation rates corre- 
sponding to random relaxation representing the macroscopic behavior of the 
material. 

From the relaxation function one can then derive the relaxation response 
f{t) = — The one-sided Fourier transform 

POO 

X{uj) = x'iuj) - tx"{to) = / exp{-iut)fit)dt (11) 

^0 
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yields the susceptibility response in the frequency range .If a = 1, the deriva- 
tive of the Laplace transform eq. (9) gives: 



/-.w^n^ifi^.i.^i (12) 

r is the Euler gamma function and U {a, b, z) is the confluent hypergeometric 
function : Jq°° e~^*t'^~^(l + t)^~"'~'^dt. Using the asymptotic behavior of 
functions f/(a, 6, z), we get for t — > oo a power law. 

lim fa=i,q{t) cc (13) 

t^oo 

when t —> 0, one gets a simple Debye behaviour. This first result indicates 
that the high t power law behavior is associated with the Tsallis q entropy 
parameter. When < a < 1, an analytical expression is not available but 
an analysis of the results for some rational values of a and implying MeijerG 
functions allows us to conjecture the following asymptotic behavior confirmed 
numerically for any couple of values (a, q) in the proper physical range. 



lim/„,,(t) oct°-^ (14) 
lim/„,,(t)ocr^-"f^ (15) 

We see that in this model, the empirical exponent n and m are simply 
related to the two parameters a and q : n = 1 — a and m = a^Ef- If a = 1, 
one recovers the asymptotic behavior of eq.l3. For intermediate time, the 
response function can be fitted to the derivative of a stretched exponential 
function with exponent a over several orders of magnitude. The condition 
< m < 1 defines the range of q : ^^j^^^ < q < 2. The narrower range 
observed in most experimental data 1 — n < m < 1 defines a narrower 
physical range: < q < 3/2. 

At this point it is of interest to compare our results with the only self- 
contained mathematical theory which can account for the relaxation uni- 
versality in the framework of modern stochastic theory fn\. It has been 
demonstrated mathematically jni^l^O] that a stretched exponential relax- 
ation function can be obtained if and only if the relaxation rate probability 
distribution is a completely asymmetric Levy distribution function. In that 



7 



case the low t exponent can be related to the Levy-a parameter i.e. n = 1— a, 
but this model does not yield the large t power law [201 • This results has 
been expressed as a consequence of the Levy central limit theorem. By gener- 
alizing further the Levy central limit theorem assuming a randomness in the 
number of effective relaxation channels , Weron et al. ^T] were able to derive 
using precise stochastic theory arguments a general simple analytical form of 
the relaxation function which reconcilliates and unifies the two power-laws 
and the stretched exponential empirical expressions. 
The Weron expression of relaxation function reads 

k) = [l + k{uj^t)"]-^/'' with 0<a<l, k>a (16) 

The asymptotic behavior of the response function f^{t,k) = —^^—iiiM fQ^- 
small and large t is given by 

{au!c{uJct)°'^^ for ujct « 1 
(17) 
auJck-^-^/''{uJct)-'^^''~^ for u^t » 1 

This corresponds with the results of our relaxation (a, q) model if we make 
the identification 

By Fourier transform it can then be shown straightforwardly from the two 
asymptotic behavior in the time domain that both models obey the Jonscher 
universal laws : 

lim = cot(n7r/2) (19) 

and 

lim , ^ ^""^ , = tan(m-) (20) 

with n = 1 — a, m = a/k, and k = (|5~) -These relations have a thermody- 
namic content pl| and the observation that m is related to a nonextensive 
entropy parameter has to be further analyzed. 
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4 Numerical results and discussion 



In Fig.l, we show the behavior of the relaxation function (f)a,q{t) for a = 0.5 
and for q varying in the range 1 < q < 2. The slowing down of the relaxation 
process increases with q. 

Fig. 2 shows the response function fa,q(t) in a log- log plot for a typical 
value ofn = 0.48 {a = 0.52) , m = 0.68 (g = 1.433). for our model, and for the 
Weron response function f^f.{t) with k = |5| iii the same time units {ujp = 
1). These values correspond to a typical experimental example quoted and 
discussed in ^U] .The two models can fit the relaxation and dielectric behavior 
of a large number of data reported in Jonsher |3] for dipole relaxation systems. 
For values around t = 1, f{t) behaves as the time derivative of a stretched 
exponential over 3 or 4 decades and again most of the experimental data 
found in the literature with a ranging from to. 0.3. to 0.8 ^H] can be accounted 
for. 

The main objective of this letter was to open a new line of thinking in this 
problem and to show that general thermodynamic principles based on exten- 
sion of the concept of entropy appears to be a natural way to understand the 
ubiquity of hierarchic self- similarity and power laws in the relaxation behav- 
ior of large classes of materials. Moreover as demonstrated in [22] the power 
law behavior of statistical distribution is rooted in the central limit theorem 
and a complete theory of relaxation based on nonextensive arguments should 
be based on formal stochastic arguments. In this letter we have shown that 
it is indeed the case since we found a simple relation between the Tsalllis 
parameter q and the stochastic parameter k in the Weron generalized relax- 
ation function. It should be emphasized that the Tsallis generalization is not 
unique. Several alternative definitions of entropy has been elaborated [ZI|H]- 
The relation between k and q is therefore not unique [Z][H1 and a general 
discussion of the link of relaxation with nonextensivity will be presented in 
further publications. We finally want to stress that the ideas developed in 
this letter are not restricted to dipolar relaxation. 
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FIGURE CAPTIONS 

Figl. Relaxation function 4>a,q{t) for values 1 < g < 2. and a = 0.5. 

Fig2. Response function fa,g{t) = -^^^ for a = 0.52, q = 1.433 (full line). 
The dotted line is the Weron f^{k, t) for a — 0.52, k — 0.765. The exponents 
n and m are.n = 48 and m — 0.68. We have used the relation k — 
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